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Abstract 

The constraints on the time variation of the fine structure constant at recombination epoch 
relative to its present value, Aa/a = {a^cc — anow)/onow; are obtained from the analysis of the 
5-year WMAP cosmic microwave background data. As a result of Markov-Chain Monte-Carlo 
analysis, it is found that, contrary to the analysis based on the previous WMAP data, the mean 
value of Aa/a = —0.0009 does not change significantly whether we use the Hubble Space Telescope 
(HST) measurement of the Hubble parameter as a prior or not. The resultant 95% confidence ranges 
of Aa/a are -0.028 < Aa/a < 0.026 with HST prior and -0.050 < Aa/a < 0.042 without HST 
prior. 



The variation of the fundamental physical constants is a longstanding issue. Dirac first 
considered such a possibility and proposeded that the Newton constant should be 

inversely proportionaol to time. While his claim is not compatible with the current observa- 
tions, recent unification theories such as superstring theories naturally predict the variation 
of the fundamental constants Sj. Because of these theoretical motivations, it is important 
to measure their possible time variation observationally. 

Among various fundamental constants, the time variation of the fine structure constant a 
has been most extensively discussed in observational contexts. We briefiy summarize those 
terrestrial and celestial limits on a as follows 

• The atomic clocks constrain the current value of the teniporal derivative of the fine 
structure constant as aja = (-3.3 ±3.0) x lO^^V"^ flflQQ. 

• a/a = ( — 1.6 ± 2.3) x 10~^^yr~^ from the measurement of the frequency ratio of 
aluminium and mercury single- ion optical clocks 

• Aa/a = (-0.8 ± 1.0) x 10"^ or Aa/a = (0.88 ± 0.07) x 10"^ from the Oklo natural 



reactor in Gabon {z ~ 0.1) 10|. 



Aa/a = (-0.57 ± 0.11) x 10"^ {z ~ 0.2-4.2) |ll|, and Aa/a = (-0.64 ± 0.36) x 



10^^ {z ~ 0.4-2.3) 12, ll3| from spectra of quasars, the former of which is from 
the Keck/HIRES instrument, and the latter from the Ultraviolet and Visual Echelle 
Spectrograph (UVES) instrument. 

-5.0 X 10"^ < Aa/a < 1.0 x 10"^ (95%C.L.) from big bang nucleosynthesis (BBN, 
z ~ lO^-lO^O) y. 

-0.048 < Aa/a < 0.032 Q, -0.06 < Aa/a < 0.01 Q or -0.039 < Aa/a < 0.010 
17( 1 (95%C.L.) from the cosmic microwave background (CMB, z ~ 10'^), the former 

two of which are based on the analysis of the 1-year WMAP data and the last one on 

the 3-year WMAP data. 



We also note that the seasonal variation effect on a has been also discussed in 

In this paper, we focus on the CMB constraint on a from 5-year WMAP data, finding 
new limits on its value at the recombination epoch. While the other physical constants may 
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vary in time simultaneously, they are so model-dependent [19| that we consider only the 
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2l| 



variation of a here. An example of such a class of models can be found in 

Both CMB and BEN 22(] are useful in order to obtain the constraints of the variation of 



a over a cosmological time scale. Although BEN is surperior in that it can probe a longer 
timespan, it has a drawback that the effect of a in Helium abundance is model dependent 
so that we cannot get a robust result from EEN analysis. On the other hand, the physics of 
the CME is much simpler and well understood with high precision data, so we can obtain a 
meaningful limit on the variation of a from the CME data. 

As is well known, changing the value of the fine structure constant affects the CME power 
spectrum mainly through the change of the epoch of recombination 23, 2^. Hence it probes 
the value of a at this particular epoch. Let us summerize the main part of the recombination 
process and see how a appears. 
Following the treatments of 
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26j, which are implemented in the RECFAST code |27l |. 



recombination process is approximated by the evolutions of three variables: the proton 
fraction Xp = rip/riH, the fraction of the singly ionized Helium XHeii = "^Heii/'^-H, and the 
matter temperature Tm- Here, nn is defined as the total Hydrogen number density. Their 
evolution equations read 



dz H{z){l + z) 

'^a;Heii _ 
dz 

dTM 



anXeXpTiu - /?h(1 - a;e) exp 



H{z){l + z) 



^bTm 

^Heia^Hell^^e'^H " PueifHe ~ Xnell) 6Xp 

2Tm 



knT, 



Xp 



M 



Tr) + 



(1) 

(2) 
(3) 



dz 3H{z){l + z)mel + fue + Xe ^ 1 + z 

In the above equations, H{z) is the Hubble expansion rate at the redshift 2, ax = 
2a'^h'^ /{STTmlc^) is the Thomson scattering cross section, aji = k'^/ {1207T c^h^) is the black- 
body constant^, z/h = c/(121.5682 nm) is the Ly-a frequency, z/He2is = c/(60.1404 nm) is the 
He2^s-l^s frequency, Xe = ne/rin = Xp + XHeii is the free electron fraction, Tr is the radiation 
temperature, and /ne = ^/(4(1 — Yp)) is the number ratio of Helium to Hydrogen, where 
Yp is the primordial Helium mass fraction which we take 0.24. Ch (Chc) is the so-called 



^ h is the Planck constant, ks is the Bohzniann constant, c is the speed of hght. 
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Peebles reduction factors 



Cu 

Cue 



1 + KnAnrinjl - Xp) 

1 + -ft'HeAHe'^H(/He 



(4) 
(5) 



1 + Kue^Aue + /5He)?^H(/He " Xnell) exp {-hUps/kBTM) 

where Ah is H 2s- Is two-photon decay rate, Ahb is Hel 2^s-l^s two-photon decay rate, 
'ps = ^'Hc2ip - i^Rc2^s, Ku = c^/iMv^^pH), and i^He = cV(87rz/|[^2ipiJ). Oh and ctHe are the 



z/, 



case B recombination coefficients and they are well fitted by 

at'' 



1 + cf^ 



Q^He 



1 + ^ 
Ti 



i+p 



n -1 



[m^'s-i] 



(6) 
(7) 



with t = Tm/IO^ [K] , a = 4.309, 6 = -0.6166, c = 0.6703, = 0.5300, F = 1.14 
and q = 10-l^■^^^ p = 0.711, Ti = 10^-"^ [K] , T2 = 3 [K]. Finally, /3h and pne are the 
photoionization coefficients 



2'KmekBTj 



M 



/5] 



He 
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3/2 
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exp 



ksTM 

^^He2ls 
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(8) 
(9) 



Now, we show how the above quantities depend on a. Since binding energies scale as a 



the frequencies u also scale as a^, and -ft'H(-^He) oc ^mpi^}ie2^p) According to j24| , 

two-photon decay rates Ah and Ahc scale as a^, and the recombination coefficients an and 
aHe scale as where we adopt £ = 0.7 following 24 1. 



As already pointed out in 
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2J], the larger value of a at the recombination epoch 



results in the higher redshift of the last scattering surface. Thus, increasing a result in 
three characteristic signatures in the angular power spectrum of the temperature anisotropy, 
namely, shift of the peaks to higher multipoles, increase of the height of the peaks due to 
the enhanced early integrated Sachs Wolfe effect, and decrease of the small-scale diffusion 
damping effect. These features can be seen in FIGH] 

We constrain the variation of a using three kinds of CMB anisotropy spectra, namely, 
angular power spectrum of temperature anisotropy, Cf^, that of E-mode polarization, Cf^, 
and cross correlation of temperature and E-mode polarization Cf^ of the 5- year WMAP 



data 



28, y, 3. 



For this purpose, we have modified the CAME code 
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32| including 



7000 



6000 



5000 



4000 



3000 
2000 



1000 




10 



100 



1000 



FIG. 1: The CMB temperature anisotropy spectra for no change of a (soHd red curve), an increase 
of a by 10% (dashed yehow curve), a decrease of a by 10% (dotted blue curve). 

the RECFAST code to calculate the theoretical anisotropy specra for different values of a 
at recombination and we performed the parameter estimation u sing Markov-Chain Monte- 
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34|. 



Carlo (MCMC) techniques implemented in the CosmoMC code 

We have run the CosmoMC code on four Markov chains. To check the convergence, we 
used the "variance of chain means" /"mean of chain variances" R statistic and adopted the 
condition R-\< 0.03 

First, we consider the modified version of the fiat ACDM model, that is, as for 
cosmological parameters we take (^b/;.^, fi/jAf/i^, i^o, '^s, t, Aa/a), where fi^/i^ is the 
normalized baryon density, i^/j^/Zi^ is the normalized cold-dark-matter density, Hq = 
lOOh [km sec~^ Mpc~^] is the Hubble constant, ng is the spectral index of the primordial 
curvature perturbation, Ag is its amplitude, and Aa/a = (orec — anow)/ttnow is the variation 
of the fine structure constant at recombination relative to its present value. We have also 
analized the standard fiat ACDM model without Aa/a and compared the other parameter 
values between these two models. 

The results obtained from MCMC calculations are given in TABLE [U and FIG|2l In 
TABLE [H we present the mean values and the 68% confidence intervals of the cosmological 
parameters and FIG. [2] shows the one-dimensional marginalized posterior distributions of the 
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FIG. 2: One-dimensional marginalized posterior distributions for the parameters of the modified 
flat ACDM model (solid red curve), and for the standard flat ACDM model (dotted blue curve). 
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FIG. 3: Two-dimensional marginalized posterior distributions for IS.a/ a-ilsh?' and Aa/a-HQ. Red 
contours are with HST prior, and blue contours are without it. For each contour pairs, the inner 
one represent the 68% bound and the outer one 95%. 
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TABLE I: MCMC results on the mean values and 68% confidence intervals of cosmological param- 
eters for the two models. 





modified flat ACDM 


standard flat ACDM 
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79 1+2.6 
'^•J^-2.6 


Aa/a 


-0.0008941°:°}^! 





parameters. From these results, it can be seen that the effect of the additional parameter 
A«/a is only to increase the errors of the other parameters and the mean values of the 
other parameters in modified fiat ACDM are practically the same as in the standard flat 
ACDM. The marginalized distributions of Hq and in FIG J2] suggest the degeneracy of 
these parameters with Aa/a. 



: Hubble Key Project 
3511. that is, we have 



Actually, in the above calcutions, we have incorporated the result o: 
of the Hubble Space Telescope (HST) on the Hubble parameter Hq 
imposed a prior that Hq is a Gaussian with the mean 72 [km sec~^ Mpc~^] and the variance 
8 [km sec^^ Mpc^^]. If we do not use the HST prior, we can only get weaker constraints 
on the parameter values because of projection degeneracy. To check the effect of the HST 
prior, we show two-dimensional marginalized distributions with and without HST prior in 
FIG|3l and one- dimensional distributions in FIGIH It is confirmed that the HST prior is 
very important in order to realistically constrain the time varition of a. 

The 95% confidence interval and the mean value of Aa/a from 5- year WMAP data with 
HST prior are 

- 0.028 < Aa/a < 0.026 and Aa/a = -0.000894, (10) 
respectively. Without the HST prior, they read 

- 0.050 < Aa/a < 0.042 and Aa/a = -0.00181, (11) 
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FIG. 4: One-dimensional marginalized posterior distributions for the results with HST prior (solid 
red curve), and without it (dotted blue curve). 



respectively. Previous results from 1-year WMAP data are —0.048 < Aa/a < 0.032 and 
—0.107 < Aa/a < 0.043, with and without HST prior, respectively [l^, so our results from 
5-year WMAP data are about 30% tighter than those from 1-year WMAP data. We also 
note that for the 1-year data the mean value of Aa/a was found to be Aa/a = —0.04 
without the HST prior although it was practically equal to with it. For the 5- year data we 
have found that the mean value remains practically intact whether we use the HST prior or 
not. This may be interpreted as an indication that the observational cosmology has made a 
step forward to the concordance at an even higher level. 

Next, we take (fis/i^, QoMh'^, Hq, Us, Ag, r, w, Aa/a), where w is the dark energy equation 
of state. In addition to the 5-year WMAP data, we use the HST and Supernova Legacy 
Survey 36|] priors here. Although the possibility of the degeneracy of w with Aa/a was 
pointed out some time ago [32], we can conclude that such a degeneracy is very weak now, 
for we have both temperature and polarization data enough to constrain w and Aa/a 
simultaneously (see FIG|5]). The 95% confidence interval and the mean value of Aa/a in 
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FIG. 5: One-dimensional marginalized posterior distributions for the parameters with variable w 
(solid red curve), and w = —1 (dotted blue curve). 

this case are 

- 0.033 < Aa/a< 0.032 and Aa/a = -0.00186 (12) 

which is shghtly larger compared with the model with the cosmo logical constant. 

In summary, in terms of the MCMC analysis using CosmoMC code, we have updated 
constraints on the time variation of the fine structure constant a based on 5-year WMAP 
data. We obtained tighter constraints compared with previous results from 1-year WMAP 
data due to the inclusion of the polarization data and the decrease of the statistical errors. 
Compared with the result based on the 3- year WMAP data [l3] , where comparison between 
the cases with and without HST prior has not been made, the resultant limit are almost of the 
same order of magnitude but the mean value of ours is closer to 0. We have verified that the 
null result is favored about the variation of a and the addition of this new parameter Aa / a 
does not essentially affect the determinations of the other standard parameters contrary to 
the case of the analysis based on the 1-year WMAP data [l^. 

We have also studied the possibility of the degeneracy between w and Aa/a, finding 
the slightly relaxed limit on the latter parameter due to the addition of w, but no drastic 
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degeneracy. Also in this case, we cannot find any evidence of time varying a in tlie 5-year 
WMAP CMB data. 



Note added 



After we have finished our analysis, we bacame aware of a paper by Scoccola, Landau, 



and Vucetich [38|, who also analyzed constraints on the time variation of a using 5-year 
WMAP data. Their main focus, however, is dependence on the details of the recombination 
scenario, which they have shown to be weak. It appears that they did not incorporate HST 
prior in contrast to our analysis. Hence our paper is complementary to theirs. 
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